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Slender vortex filament with slowly varying core structure 

ABSTRACT We give a brief review of the asymptotic theory of slender vortex filaments with emphases 

on the choices of scalings characterizing the physical problems and the corresponding assumptions and/or 
restrictions introduced in the formation of the asymptotic theory of Callegari and Ting (1978) and its 
extension by Klein and Ting (1992). In particular, the slender filaments considered are assumed to be 
forming loops or tori. Because of this restriction, the theory is not applicable to the trailing vortex system 
of a rotorcraft. We describe the multiple length scales characterizing the vortex system, formulate the 
expansion scheme, derive the governing equations and then identify the assumptions or restrictions inherent 
in the multi-scale analysis and needed for the validity of the asymptotic theory of the trailing vortex system. 


1. Introduction 

It is well known that the velocity of a filament depends on the vorticity distribution in the core and 
the velocity is undefined when the filament is modeled by a vortex line of zero core radius [1]. Consider a 
slender filament, with centerline C : x = X(£, s), characterized by a parameter s. Here we choose s to be the 
arc length at t = 0. The velocity Q induced by a vortex line C is given by the Biot-Savart formula 

Q(f, x) = (T/(47 t)] J [X' - x] x dX'/\X' - x| 3 . (1.1) 

For a point x in the neighborhood C, i. e., in the normal plane of point X(t, s) on C, we write x = X + rr, 
where r — ncos<f> + 6sin0, r and <f> are the polar coordinates and r,n and b are the unit tangent, normal 
and binormal vectors of C. In approaching the vortex line C ) i. e., r = |x — X| — > 0, the formula gives three 
singular terms, [2] - [6], 

Q = \T/(2nr)\8 + [T/c/ (47 t)] ln[S/r]6 + \Tk/{±tt)\ 9 cos<f> + , (1.2) 

where Q' denotes the remainder of the Biot-Savart integral, and 6 = f x f denotes the unit circumferential 
vector. The first term with l/r singularity corresponds to the circumferential velocity of a 2-D vortex point 
in the normal plane. The second term, the binormal velocity with lnr singularity, represents the curvature 
effect. The third term, a circumferential velocity depending on 0, does not have a limit as r — ► 0. After 
subtracting these three singular terms from Q, the remainder Qf has a limit, known as the finite part of the 
integral. These singular behaviors are not valid for a real fluid because with viscosity the flow field has to 
have a continuous velocity gradient. 

The B-S formula (1.1) for the vortex line C is valid in the “outer region”, i. e., away from the filament, 
of a normal length scale t. Near the filament, of a typical core size <5, the singular terms on the centerline C 
in (1.2) shall be matched with or removed by the inner solution for the core structure using the methoc^of 
matched asymptotics. The asymptotic analyses with the small expansion parameter, slenderness ratio, 

c = 6/e, (1.3) 

were carried out in a sequence of papers [2, 3, 4, 5] from the 2-D cases in 1965 [2] to filaments in 3-D space 
with large swirling and axial flow in the core in 1978 by Callegari and Ting (C-T), [5]. The analyses were 
summarized by Ting and Klein in 1990. [6]. 

In the analyses of C-T for a slender filament we assume that: 

(i) The filament forms a slender torus, of length S(t) = 0(£) } so that all the physical entities in the inner 
region around the centerline C and the geometry of C, are periodic function of its arc length s with 
period 5, or its initial arc length s with period So, i. e., 

X(t, s 4* So) = X(£, s), and /(£,r, 6, s + So) = f(t , r, 6 , s) , (1.4) 

where (r,0, s) denote the curvilinear coordinates with respect to C and / denotes an entity near C, i. e, 
the velocity components and the pressure. 
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(ii) There i, only one Wed length scale l other than the typical core size J. This Implies that the 
orderT ltUre ’ ^ ^ ° f the dement' end the length scale of the background flow A, are 


radius 
of the 


1 /* = 0(f), S = 0{i) and A = 0(£). 


Here, the distance between two adjacent filaments is a length scale for the background 
ese two typical length scales define the slenderness ratio, 


flow to be 


(1.5) 

0(A). 


e = s/e. 


( 1 . 6 ) 


(m) There is one typical velocity scale U for the background flow and the velocity of a filament centerline is 
or the order of (7, i. e., 


X(«, s) ■ h = 0(U) and X(t, s) ■ b = 0(U) with X-f = 0, 
where f, n and b denote the unit tangential, principal normal and binormal vectors of C. 


(1.7) 


(iv) The core structure does not vary along the centerline of the filament, i. e., independent of the social 
variable s. 

Rpvrfr!^ 6 a r f flU D iS V ir7 US ’ probIem is characterized by the physical parameter, the background 
y Ids number, R e - Ut/v, where v denotes the kinematic viscosity.f We consider R e » 1 so that 
the viscous effect m the background flow is of higher order while in the core structure the viscous effect is 
important when the core size is very small initially and/or the effect is accumulated over a long time We 
choose the “distinguished limit” , 


1 /s/R e =0{t) or l/y/R e <Ke, 


( 1 . 8 ) 


so that the equations for the evolution of the core structure retain the nonlinear convection (inviscid stretch- 

ing) and viscous diffusion terms. The constant K represents a typical ratio of viscous to inviscid effect The 
inviscid limit is K — + 0. 

In 1992, Klein and Ting, [7], extended C-T theory to allow for axial core structure variation, thus 
removing assumption (iv) while keeping (i), (ii) an d (iii). They arrived at a system of Integra-differential 
equations much more complex than that of C-T. 

Because of assumptions (i) and (ii), the analyses of C-T and K-T are not applicable to the vortex 
system trailing rotor blades, for which the filaments do not form tori with length of the order of the radius 
ol curvature. The filaments have at least three distinct length scales in addition to the core size 5 They 
are the typical distance between two filaments, d, the radius of curvature f, and the length scale L for the 
axial variation of core structure and size. It is clear that the scale d is of the order of the pitch of a filament 
divided by the number of blades. The slender filaments reach a typical length L when they merge i e 
when the core size becomes comparable to the distance between two filaments. We have 


T » f, d » 5 while d « i or d = 0(f), 
We consider d = O(f) and hence has only one more small length ratio 


(1.9) 


i / L = ( 1 . 10 ) 

in addition to the slenderness ratio, e in (1.6). 

. co To 7 XP f Iain or w here the existing theory is not applicable to the rotorcraft vortex system, we review 
10 formulatlon of C - T theor y [5] based on assumptions (i) to (iv). We show the expansion schemes 

and the two systems of equation s defining the motion of the filament centerline and the evolution of the core 

t When the vortical core is turbulent, * shall be replaced by an eddy viscosity v d , resulting in a relatively 
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structures respectively and then the coupling of these two systems and the 
(ii) for the closure of the second system. 


necessity of assumption (i) and 


In j* 3, we mtro Juce the expansion schemes for filaments having three length scales S, l and L, removing 
assumption (1), modifying (11) and (in), while keeping (iv). We describe the multi-scale analysis yielding two 

new systems of governing equations and then discuss the conditions on the filament centerlines needed for 
the multi-scale analysis. 


2. Filament in the form of a slender torus of length 0(t). 

We re view the C-T theory [5] for slender filaments submerged in an incompressible irrotational flow The 
background flow in absence of the filaments is a potential flow, fulfilling the Navier Stokes (N-S) Equations 
We consider the strength of a filament, the circulation T, to be 0(U£). Let t and U be the unit length and 
velocity, i. e., I = 1 and U = 1, and hence tjU = 1 is the unit time. With the core size 5 = 0(e) we have 
the order of magnitude of the vorticity Q and that of the velocity v in the vortical core- 


|fi| = 0(e~ 2 ) and |v| = 0(e~ 1 ) . 


( 2 . 1 ) 


To construct the inner solution, i. e., the core structure near a point X on the centerline, we introduce the 
intrinsic coordinates, r, 4>, s and then replace the polar angle <j) in the normal plane by the angle 


0 = <j>- 9 0 (t, s ) , with dOo/ds = -aT, a = |X,|, 
so that the new coordinates, r, 9, s, are orthogonal with 

— f dr + 9 rd9 + t h^ds where /13 = <rX, A = 1 — rn cos <}> . 

The rotation 9o(t,s) from <fr to 9 accounts for the torsion T(t,s) of C. See [4]. 

For the inner region, i. e., r = 0(e), we introduce the stretched radial variable 


( 2 . 2 ) 


(2.3) 


t ~r/e and assume X t = 0(1) with X t • f = 0. 

With the assumption (iii), we introduce the velocity V(t,r,9,s) relative to C and write 


(2.4) 


v “ X ‘ + V w ith, V = ur + v9 + wt = 0(e x ) and u = u = 0 at r = 0. (2.5) 

The continuity and momentum equations in the moving curvilinear coordinates are 

(rAu) r + (Xv)e + (w a + X st ■ t)(r/cr) = 0, (2.6) 

dV/di + [- ~ r7VT ]X t + X* = -Vp + t'AV + ^[^X s ], (2.7) 

To characterize a slender vortex filament, or to be consistent with (2.1), we need a large circumferential 
velocity, v = 0(e ), and allow for a large axial flow w, but the radial component u remains 0(1) Note 

that the inner solutions are periodic in 9 with period 2 tt and in s with period S 0 on account of (i). Using 
this physical model, we arrive at the expansion scheme, 


f(t,r,9,s,e) — e m {f (0) (t, r, 9, s) + e/ (1) (t,f, 9,s) + 0(e 2 )} . (2.8) 

We consider ln(l/e) to be 0(1) with respect to e, therefore, f^\j = 0, • • • , could be functions of ln(l/e). If 
/ stands for v, we have m= 1 , with u °) = 0. If / stands for p, we obtain m = 2 from the radial momentum 
equation. We substitute the power series (2.8) into the N-S equations, equate the coefficients of like powers 
of e, and obtain systematically the sets of leading and higher order equations. 
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With u<°> = o, we obtain from the N-S equations, (2.6) and (2.7), the leading order continuity equation 
and axial and radial momentum equations, 

Vg 0) = 0, w ( e 0) = 0 and p ( f 0) = [u (0) ] 2 /f . (2.9) 

Then the circumferential momentum equation becomes an identity. The solution, independent of 0, represents 
a quasi steady two-dimensional axi-symmetric flow in the normal plane of C. 

Note that the dependence of the core structure, v (0) and w (0) , on t,f, and s is yet to be defined later 
by the compatibility conditions of the second order equations (5, 6, 7]. 

The next order (the first order) system of equations is 


v + (ru (1) ) f + rw {0) I a + rv W K sin (f> = 0, 

+ Wff^v^/r + w^v^k sin</> = — (p^ + w^w^^/a, 
+ « (1) (fw (0) )f - (w (0 ^) 2 f«sin</> = -p ( / } - fw(%(V/a, 
V (0) Ug^ - 2t/°V 1 ) + (w^ 0) ) 2 fK COS0 = -fp^K 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


They are linear equations for the first order solutions, u^\v^ etc, with inhomogeneous terms nonlinear 
in the leading order solutions. We can decompose the first order solutions and their governing equations to 
their axi-symmetric and asymmetric parts in 0. For example, 


f(9) =< / > +/ o (0) with 


< / >= ^ ' / d9f(6) and < f a (9) >= 0 . 


(2.14) 


Here < / > denotes the 0-average of / and f a denotes the asymmetric part of /. Likewise, we can do the 
decompositions for the higher order solutions and their governing equations. Because of the distinguished 
limit (1.8) and our assumption of only one time scale i/U, the viscous terms and the time derivatives do not 
appear in the first order equations (but will appear in the second order equations). 

The 0-averages of those first order equations yield a system of quasi-steady inviscid equations for the axi- 
symmetric parts of the solutions. Consequently, the core structure for t > 0 has to fulfill two be consistency 
conditions identified with two classical relationships for steady inviscid flows. See [7], Early on, C-T [5] 
found a special case fulfilling both relationships. It is: 


v™ = = 0 an d p£ 0) =0 if any one of them holds. 


The C-T theory is based on the restriction that 


= 0, hence = 0 and p£ 0) = 0. (215) 

Thus the leading order solutions, and p (0 '> are functions of t and f only. 

By matching the asymmetric inner solution of the core structure asf-»oo with the outer solution for 
T — * °> tfie second and the third singular terms in the Biot-Savart integral (1.2) are removed and the velocity 
of the centerline, X t , is defined by 

X t (f, s) = n[Q 0 • n] + 5[Q 0 • b] + S[r«/(47r)][ln(l/e) + C v + C w ] , (2.16) 

where Qo(f,X) = 4- denotes the background velocity without the filament plus the finite part of Q 

in (1.2). C v (t,s ) and C w (t,s ) denote respectively the contributions of the circumferential and axial velocity 
in the core to X £ . They are, 


Cv = f !™ (^T J Q f ’ v2 df/ ~ ln f ) + \ - c *> - J q r'w 2 dr' . (2.17) 

With ln(l/e) considered 0(1), we treat the terms in (2.15) as 0(1) while terms O(e) have been omitted. 
Their omission implies that the superscripts (0) for X and its geometrical entities, o, k, T, t etc have been 
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suppressed, i. e., we are considering X(t, s,e = 0). The terms omitted can come from the higher order core 
structure and from the effect of nonzero core size not accounted for in the B-S formula (1.1). 

The asymptotic theory of C-T requires that the initial core structure should fulfill the consistency 
conditions (2.15), and the initial velocity of the filament centerline, if assigned, has to agree with (2.16). 
These restrictions on the initial data are the results of having only one time scale. In case the initial data 
violate any one of those restrictions, we need multi-time analysis. See [2], [4], [61. 

The second order equations involve the first and second order unknowns, z/ 1 *, z/ 2 \ etc. To remove those 
unknowns, we first remove the second order unknowns, and . by using the periodicity of the solution 
in the variable 8 with period 27 r. Carrying out the 0-average of the second order continuity equation and the 
axial and circumferential components of the second order momentum equations, we arrive at the following 
equations: 


i(f < ,P'< >), + ~\< W<‘> >, +j<“>], or 


< >- r < >. 


f'df'- 


w, 


2a<°>’ 


(0) ~ 7-( f 4 0) )r = - Jo)( w(<>) < WW >* + < p(1> >•) 

(j(0) w r <u > 


(2.18) 


(2.19) 


and 


lf (0) 


- Kr 


l(r4 0) )f - ^ (0) ] =- 


W 


( 0 ) 


< v (1) >, 


(7(0) 

- (rt) ( °)) f < «W > /f 


( 2 . 20 ) 


Here K 2 = v/(Uit 2 ) = 0(1) on account of (1.8). 

To remove the first order unknowns in (2.18), (2.19) and (2.20), which are the s-derivatives of their 
^-averages, we make use of assumption (i), the periodicity of s with period S 0 expressed by (1.4). We use 
(2.18) to remove < u ( 2 ) > and then integrate (2.19) and (2.20) with respect to s over the entile length S(t) 
or s over Sq with ds = c^ds. We arrive at the compatibility conditions of the second order equations, 
which in turn are the evolution equations of the core structure, v^(t,r) and w^(t,f). They are: 


wl 0) = (tf 2 /f)[f4 0) ] f + (f 3 /2)(5/S)[u/°Vf 2 ], , 



With the axial vortipity C (0) related to n (0) by fC (0) = (rv^) f , (2.22) is then replaced by 


( 2 . 21 ) 

( 2 . 22 ) 


C f (0) = (K 2 /f)[r(j.% + (5/ S) (r 2 £(°))f / (2r) . (2.23) 

We use the axial vorticity as a primary variable instead of the circumferential velocity v (0) , because C (0) 
decays exponentially in f, while becomes r/( 27 rf) as r — ► oo. Thus the inner solution vo matches or 
removes the leading singular term of the B-S integral (1.2). Recall that the second and third singular terms 
were removed by the asymmetric part of the first order solutions, therefore, we have shown that the leading 
and first order inner solutions have matched or removed all three singular terms of the outer solution (1.2). 

Equations (2.21) and (2.23) can be considered as the evolution equations for w^(t y s) and C^(^> 5 ) 
respectively, but with coefficients depending on the filament length S(t) ) they are coupled with the equation 
of motion (2.17) of the centerline, X(£,s) , and hence with each other. The radial integrals of the core 
structure, C v and C w in (2.17) in turn couple the equation of motion with the two evolution equations, of 
the core structure. Together, (2.17), (2.21) and (2.23), form a close system of equations for the dynamics 
of filaments fulfilling assumption (i), i. e., in the forms of slender torii. Solutions of this system for a given 
set of initial data were presented in [5] and [6]. In the next section, we shall explain how to derive a close 
system of equation for slender filaments without assumption (i). 
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3. Filament trailing a rotor blade - Filament with multiple length 
small core size 5 ♦ 


scales in addition to the 


The vortex system trailing a helicopter rotor is extremely complex. It can be treated as a steady axi 

TdSrof * for r d the fiow is unste8dy “ d “is 

g the flow field wlth rotor tip speed nearly transsonic and large forward velocity is bevond reach 
m the near future. Here we shall describe how to modify the theory of C-T for flows whh sleX v^x 
filament [5, 6] to simulate the vortex system trailing a rotor. Note that the theory of C-T allows for unsteady 
asymmetric flow but incompressible. The last restriction implies that the flow field is at low Mach number 
s in icated in § 1 , we need to consider the length of a filament much larger than the reference radius of 
curvature l of the filament and the filament structure varying along its centerline in a length scale L >> 7 
introduced in 1.10) The assumption (i) for the C-T theory shall be removed and be replaced by that for a 
two length scale analysis. See for example [8] and [9], * 

3.1. Multiple length scales and the expansion scheme 

We replace the axial variable s by two scaled variables 


s = s/e and £ = s/L . 


(3.1) 


A function / (s) is expressed as a function of the two independent scaled variables, F(s, f) with the 
tion of the existence of the average of F in S over a large interval A, but A <1 /c, i. e.’, 


assump- 


f(s) = F(s,f), when £ = ss, and 

ft A+A/2 ,£ /<+A/2 


rr r A + A/2 , /-fA+A/2 

'™ ll L '«)*■')/[/ . 

•/fA-A/2 Ji/f-X / 2 


(3.2) 


With a - ds/ds Where s denotes the arc length of C at instant t, the average of / is the line integral along 
a segment of C divided by its length, which is much longer than l but much shorter than L. Note that the 

independent variables t,r and 9 in /, F and F have been suppressed. From (3.1) and (3.2), we have 


^ d 9 f — d$F -j- and a l d s F = d$F = 0 . ( 3 . 3 ) 

For a filament / stands for the core structure variables, u, w, <, and the centerline position vector X. Instead 
of introducing new symbols for the corresponding functions of two scaled axial variables, we consider v w ( 
and X as functions of the two axial variables. This is equivalent to replace F by f and s by s noting £ = l. 
The expansion scheme (2.8) for the inner solution in one axial variable is replaced by: 


u{t, f , 9,s,£,e) — u^(t, r, 9 , s, £) + + • • • , 

v(t, f , 9, s, £, e) = e-V 0) (t, f , £) + v (1 > (t, f , 0, s, f ) + ev (2) + ■■■, 
w(t, f , 9, s, £, e) = e- V°>(f , f , 0 + w<« (t, f , 0, s, 0 + ewM +■■■ 
p(t,f,9 lS ,t,e) = e~ 2 pW(t,r,0 + e- 1 p (1 \t,f,9,s,Z)+pW + . .. 
C(t, f , 9, s, £, e) = c" 2 C (0) (t, f , 0 + e- 1 ^ 15 (t, f , 9, s, f ) + C (2) + ■ • • 


(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 


X(f, s, 0 = s, 0 + eX<« + • • ■ . (3.9) 

<t(«,s,0 = X, -T = a (0) (t,s,O . ( 3 . 10 ) 

For the core structure, the leading order solutions have only the slow axial variation in £, but both axial 
variables, s and £, can appear in the higher order solutions. For the velocity of the centerline, the leading 
term X 0) has to depend on both s and f because its radius of curvature is 0(e). From (3.2), we know 
t at a £- derivative term will not appear in the system of equations where s-derivative terms first appear. 
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Thus the leading order and first order systems of equations, (2.9) to (2.17), in §2 remain valid with the long 
axial variable £ appearing as a parameter. The removal of assumption (i), requires special attention to the 
evaluation of the finite part Q / of the B-S integral ( 1 ) at point X(s) for s > 0 with t treated as a parameter. 
We need to assume that 

(v) Far downstream, a trailing vortex centerline should be of the shape, e. g., a circular helix, which renders 
the B-S integral finite as its upper limit tends to infinity. 

(vi) Let X( 0 ) denote the starting point of the tip vortex, trailing from a rotor blade near its tip at distance 
b from the axis of the rotor. Then we need to model the vortex line for s < 0 in order to evaluate the 
B-S integral. A simple model is having the vortex line bounded or rotating with the blade for $ < 0 
until the root of the blade at distance r 0 from the rotor axis, where s = —(6 - r 0 ) and continue for 
s < —(6 — ro) along a free “root” vortex line, trailing from the root to downstream as s — * — oo. 

This root vortex line has to be present so that the circulation vanishes along any a contour around the 
trailing vortex system. This mathematical model could be realized in a experimental setup for a single rotor 
blade. In a real problem, the root vortex lines would interact with the body. In the mathematical model, the 
radius of curvature of the “root” vortex is O(r 0 ). For a small root radius ro <6 = 0(1 ), the contribution 
of the small radius of curvature of the root vortex to the B-S integral at a point near the tip vortex can be 
ignored. Thus we can model the “root” vortex by a curve with radius of curvature o{t). In a hover flight, 
the root vortex line can be approximate by a straight line, along the rotor axis. In this special case, the 
vortex line for s < 0 is composed of a finite line segment and a semi-infinite one, and explicit formula for the 
B-S integral for s < 0 is available. See for example [ 10 ]. 

3*2. System of governing equations for “long” filaments 

To study the evolution of the core structure in the long length scale L, we introduce the distinguished 
limit, 


C = 0(e) , i. e., the constant G = <;/e = 0(1) . 


(3.11) 


Then a Gvf term shall appear in the second order equations together with a vi X) term. The same applies 
to w and p. With the addition of the ^-derivative terms to (2.18), (2.19) and ( 2 . 20 ) in § 2 , they become 
respectively, 


i(f < u™ >) f + i[< u> (1) >* +Gw { V + <x (0) ], or 


< u< 2 > > 


L [< ™ (1) >s +Gw ? ) ] f ' df ' 


f«r<°> 
2<7<°> ’ 


w t 0) ~ ^r(r4 0) )f = ^r(r < u (2) > f --L[< p W >s +Gp< 0) ] - 


„(°) 




< u ^ 


and 


v { t 0) -K> 


J(r4 0) )f - ^ (0) 

Note that from (2.9) we have 


= ~~(oj'[ < >* +G^ 0) 1 - (rv {0) )f < u (2) > /f. 


( 0 ) 


rOO 


If ■ 


(3.12) 

(3.13) 

(3.14) 

(3.15) 


To remove the first order unknowns < t^ 1 ) > 3 and < vv > s in (3.12), (3.13) and (3.14), which are the 
s-derivatives of their ^-averages, we carry out their s-average in the sense of ( 3 . 2 ) while making use of ( 3 . 3 ), 

^ s 

we arrive at the equations for < u (2) >, w (0) and v (0) which are functions of t,r and £. They are: 


-(r < U W >) f 4 . Gwi + 


5 (°) 

5W 


or < u 


( 2 ) 


- ji : 


4 0) + Gw ( 0 ) w^ 0) + < u (2) > 4°) = {K 2 /f)[rwf ) ] 


it^Vdr' - tzl, 

(3.16) 

— rw^E , 

(3.17) 
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and 


where 


v t 0> + Gw (0) t)| 0) +^/^(ft,(°>) r -/f - AT 2 [i(f4 0) )f - 4u (0) ] , 

H(f , 0 = = Bm {[ £ /t+V2 <r t (0 >(<, f, £)<&']/[ / ?A+A/2 o-f°>(t, 5, 0*1} 

°° ^C/t-A/2 Ji/t-X/2 


(3.18) 

(3.19) 


denotes the average of the rate of change of the arc length C in the sence of (3.2) with F(t,s, 0 = 

a t (^> 0/°^ ^ ■ Using (3.16) and (3.15) to eliminate < xS 1 ' > > and respectively, (3.17) and (3.18) 
become integro-differential equations for u (0) and w<°>, in t > 0, f > 0 and £ > 0. The unknown ~(t £) 
provides the coupling with the velocity X<°>(t,s) of C, given by (2.16) in variables, t > 0, s > 0, with the 
variable £ in C v (t.£) and C w (t,{) identified as sq. 


Specify the initiaI data for the figment centerline X(0,s) and the core structure 
v ^ and w Note that the initial centerline of a filament simulating a rotor blade has 

to fulfill the assumption (v) far downstream. The centerline for s < 0, has to fulfill asumption (vi) for t > 0 
X U’ S ) = Kfcoswt + Jsinut) for r = -s e [r 0 ,6]. In addition, the initial centerline has to 
tu hll the basic assumption that the average define by (3.19) exists for A » 1 but < !/<;. We then have to 
solve numerically the system of equations and find an upper estimate of time t c > 0, below which the basic 
assumption of two-scale analysis (3.2) remains valid. 


3.3. Modeling initial filament centerlines 


Finally, we shall present a few admissible initial centerlines in the form of helices. See for example [11], 
Let k represent the axis of a helix and a e [0, 7r/2] denote the constant angle between the axis and a tangent 
line of the helix. Here we have k pointing downward. Then its principal normal vector lies in the xy plane 
and k/T = tana. We can express the unit principal normal vector n and its orthogonal unit vector rh in 
the xy plane in term of the angle <p as 


h — — sin 4>i -{- cos j , m — cos <f> i 4- sin <f> j , 
and then the unit tangent and binormal vectors of the helix, 

f = cos q H sin a m , b = sin a k -f cos a rh . 


With rh' = n, we have 

= sin ah = s'(0)/cn hence a = s' = sinaR(<f>). 
where R = 1/k. ’biven R(4>), we have the torsion T{d>) = cot a/R{cj>), the arc length, 


s{(j>) = sin a [ R(4>')d<p' , 

Jo 


(3.01) 


(3.21) 

(3.22) 


(3.23) 


and the centerline 


r<t> 

X(^) - x(0) = sina / d<p' R{(j>')[cosak + sina(cos^'i + sinfl/j)] . 

J o 


(3.24) 


FYom (3.24) we see that X(^) is a linear function R(<f>). We can get different helices by linear combina- 
tions of different R{<j>)' s. 

§3.3.1. Circular Helices. Let R(4>) = R x = constant. We could identify Ho as 6. The centerline for 
s = (osina)^ > 0 is 


Xi(*) = Hi sina[fc0 + sina(£cos<£ + jsin<£)]. 


(3.25) 
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We 3 hL^(o”" g fif+ r c U i5 SET. >‘ot C< “ P< Wh " e ° “ d " ^ P ° Si “ Ve C °" 5,antS ' 

X 2 («?!i) =Xi( 0) + cfcsina[l - e~ q ' t ‘)/q + 

[c cos Qf/(1 + 9 2 )]{i[q + e -, *(sin0 — q cos </>)] + j[q - e~ q<t> [cosf> + gsin0]} . ^ 3 ' 26 ^ 

§3.3 S. Helices in forward Motion without Tapering. Consider helix in forward motion with velocity 
Ui, being circular when U = 0. Let R 3 (J,) = R 1+gC os<t>. We have for <f> > 0, Y 

X 3 (» = Xi (4>) +g{k sin q sin ^ + cos a[i(2<f> + sin 2f>) + j(l - cos 20)]/4. ( 3 . 27 ) 

It IT' 63 ^ ^ With Ped0d P = 27r/w then X ' * advances ^ 2 *9 = UP. Thus we set g = U/u. It is 

clear that we can create a tapering helix in forward motion and other types of helices. 
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